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Abstract 

Several methods to improve the parton-shower description of hard processes 
by an injection of matrix-element-based information have been presented over 
the years. In this article we study (re) weighting schemes for the first/hardest 
emission. One objective is to provide a consistent matching of the POWHEG 
next-to-leading order generator to the Pythia shower algorithms. Another is 
to correct the default behaviour of these showers at large transverse momenta, 
based on a comparison with real-emission matrix elements. 
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1 Introduction 



With the start of the LHC in mind, there has been a recent focus on improving the descrip- 
tion of event topologies and cross sections, going beyond the Born level for many processes 
of interest. Firstly, it involves the effects of events with one or more extra jets in the final 
state, which may affect the impact of background processes and thereby the choice of anal- 
ysis strategies. Secondly, it includes next-to-leading order (NLO) corrections to production 
cross sections, which are needed for precision tests of the Standard Model and, hopefully, 
of physics beyond it. 

For the first point, real-emission matrix element (ME) calculations give a good descrip- 
tion of hard and widely separated jets, while parton shower (PS) models give the correct 
behaviour in the soft and collinear regions of phase space. The goal is to find a way to 
combine these two methods, so that each is used in its region of validity, with a smooth 
transition between the two in all physical distributions. This is not a trivial task. One 
key issue is that ME calculations describe inclusive events, while the PS generates exclu- 
sive ones. The CKKW method solves this by using ME's supplemented by analytical 
Sudakov form factors to go from an inclusive to an exclusive language in the hard region, 
and then switching to a PS below some ME cutoff scale. In CKKW-L [2] the Sudakovs 
are instead generated from fictitious showers using the same PS algorithm as in the soft 
region, to improve the consistency and continuity. These methods continue to evolve [2111] 
while other approaches include MLM [S] and Pseudo-showers [6J . Comparisons between the 
methods have been made [ZIIH]- 

For the second point, the virtual correction terms required at NLO make the calculations 
more complicated. Analytically the cancellation of real and virtual ME divergences occur in 
the soft /collinear region, i.e. where we would rather want to use the PS description. The first 
approach to solve this issue for nontrivial cases was MC@NLO [Hllin]- Here the analytical 
expression is derived for the phase space population by the first shower emission, in the 
absence of a Sudakov form-factor correction. The difference between the real-emission ME 
and this analytical PS expression, which should be finite in the soft and collinear limits, 
defines the differential cross section for events with one real "ME-based" emission, from 
which the shower should start. The rest of the cross section, wherein the analytical PS 
divergences and the virtual divergences cancel, gives the events where the shower should 
start from the lowest-order process. 

The MC@NLO approach has the disadvantage that the PS emission rate may well be 
higher than the ME one in some parts of phase space, in which case one is forced to introduce 
negative-weight events. Furthermore, the analytical shower expressions are specific to a 
particular PS algorithm, so lengthy work has to be redone not only between different 
generators but also for minor changes inside a given generator [TT] . Both of these issues are 
solved with the POWHEG method p^[T3] . where the ME's themselves are exponentiated 
to provide a pro cess- dependent Sudakov form factor. Thereby, a positive-weight algorithm 
can be obtained wherein POWHEG (almost) always generates one emission, chosen to be 
the one with largest transverse momentum. It is then up to the subsequent shower to 
respect this constraint, but otherwise without the need for a tight connection between the 
ME and PS stages. 

When POWHEG is used with showers that are not p_L-ordered, then further thought 
must be given to the interface. Specifically, in HERWIG, with its angular ordered showers 
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[imS], the first emission is at largest angle but not necessarily at largest p^. This leads to 
the idea of a "truncated shower" [12], where the shower is modified such that the hardest 
emission is generated with a modified Sudakov form factor. Subsequent emissions are then 
generated with the usual algorithm, but with a veto so that they cannot be harder than 
the hardest emission. 

The POWHEG approach is especially convenient, then, if the shower itself is p^-ordered, 
as is now the case in the Pythia generator p^HTS] . Nevertheless there are some subtleties 
that should be taken into account to optimise the interface. We discuss these issues in 
Section 2, and in particular interface to POWHEG-hvq flQ], an event generator for heavy 
quark production in hadronic collisions at NLO QCD, where some simple comparisons are 
presented for top and bottom production. We also introduce a "poor man's POWHEG" for 
cases where NLO calculations are available, but only in the traditional phase space slicing 
approach, where shower Sudakovs can be used to provide a smoother matching. 

One should note that the more sophisticated the description aimed for, the higher the 
price in terms of manpower that goes into the detailed simulation of a specific process. 
The point of injecting ME information is precisely to move away from the universal shower 
behavior, which means that each new process must be considered "from scratch". To 
provide a sensible first estimate, however, it is useful if the shower can be improved to 
get at least the qualitative behaviour right for a wide range of processes. In Section 3 we 
use the MadGraph/MadEvent [20] generator to look at further pair-production processes 
with a jet in the final state. By adding a Sudakov to the cross-section for real emission 
of a jet, to approximate the prescription used by the POWHEG method, we are able to 
make comparisons to the Pythia shower, and find modest changes that improve agreement 
significantly. 

Specifically, we will address the issue of "power showers" vs. "wimpy showers" . In |2T] , 
the authors compare two extreme choices for the maximum emission scale of the parton 
shower; either the full CM energy of the incoming hadrons (power) or the transverse mass 
of the particles produced in the hard collision (wimpy), for both virtuality- and p_L-ordered 
showers. Their conclusion was that these options bracketed the matrix-element behaviour 
for the top and SUSY production processes studied, but also that the spread of "predictions" 
is large. 

One must note that an ultimate goal would be to have a matching scheme that allows 
for both a matching to multiple real emissions and to virtual corrections. Some algorithms 
have recently been proposed in this direction [22l - [2^ . but are not yet at a stage to be used 
for serious LHC studies. We will not discuss such issues further here. 

The outline of the article is to study the POWHEG approach and its relation to Pythia 
in Section 2, to use MadEvent to gain an improved understanding of sensible default shower 
behaviour in Section 3, and to provide a summary and outlook in Section 4. 
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2 POWHEG 



2.1 The Pythia merging strategy 
2.1.1 Evolution equations 

The Pythia parton shower orders final-state radiation (FSR) emissions in terms of an evo- 
lution variable Q^, such as (previous Pythia versions) or p\ (Pythia nowadays), with 
an additional energy-sharing variable z in the branching. For QCD emissions, introducing 
t = In (Q'^Z-^qcd)^ the DGLAP evolution equations lead to the probability for a splitting 
a ^ be 

(y 

(^V = Y.^Pa^Uz)^tdz , (1) 

b,c 

where Pa-j-bc are the DGLAP splitting kernels. This inclusive quantity can be turned into 
an exclusive one by requiring that, for the first ("hardest") emission, no emissions can have 
occurred at a larger Q^. The probability that a branching occurs at t is now given by 



dP 



Xa^bcit) exp ^ - dt' J2 la^bcit') \ , (2) 



with 



Ia^bc{t) = dz— Pa^bciz) 

^ Zm.i.n. (t) ^7T 



'Zmax(t) ^ as{t) 

The introduction of this Sudakov form factor turns the unnormalised distribution into a 
normalised one, i.e. with unit integral over the full phase space. In practice, a lower cutoff, 
to, is introduced to keep the shower away from soft /collinear regions, which leads to a 
fraction of events with no emissions inside the allowed region. 

For initial-state radiation (ISR), the evolution is performed using backwards evolution 
[25] , where a given parton b entering a hard scattering is unresolved into a parton a which 
preceded it. Here, the parton distribution functions, refiecting the contents of the incoming 
hadron, must be taken into account. Such a change leads to a Sudakov with the form 

Sb{x,tmax,t) = exp - / dt }^ dz— Pa^bc{z) , , ■ (3) 

\ Jt 27r xfb{x,t') J 

One feature of the above equations is the running renormalisation and factorisation 
shower scales, i.e. the scales at which ag and the PDF's are evaluated. For both ISR and 
FSR, as is evaluated at the p± scale of the emission (the definition of p± in this context 
is discussed later in Section 12. 3p . For ISR, the fiavour dependent ratio of PDF's given in 
eq. ([3]) is evaluated at the selected t value which, for the current versions of Pythia, is p"^. 
Thus, the renormalisation and factorisation scales are the same. 



2.1.2 The merging strategy 

Probably the first use of an explicit matching of PS to ME, the so-called merging strategy, 
was introduced to handle the case of three-jet events in e"'"e~ annihilation |26]. An outline 
is given below, starting from the Born cross section for the lowest-order process e"'"e~ — >■ 
7*/Z0 ^ qq. 
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On the ME side, the real-emission cross section e"'"e — )■ ^* jZ^ — )■ qqg is given, for 
massless quarks, by the well-known expression 



ctb dxidx2 27r 3 (1 — — X'l) ' 

with cts typically evaluated at a scale of s, the invariant mass of the system. 

The DGLAP inclusive q — )• qg emission probability in (Q^, z) space can be mapped onto 
the (xi, X2) space 

^ «s(pi) 4 1 1 + d(Q^2:) 
^ 27r 3 g2 i-z d(xi,X2) ■ ^ ' 

The sum of emissions off the q and q gives W^'^ = W^^ + W^^. With the addition of a 
Sudakov form factor, as above, this becomes 

WlL^UQ') = W^'m exp J^;- dQ- W^^Q'^)^ . (6) 

For ease of notation we have here omitted the dependence of W^^ on z and the need of an 
integration dz' over a range ^min(<5'^) < z' < z'ms.xiQ'^) ^^e exponent. 

It now so happens that the Pythia shower algorithm covers the full three-jet phase 
space and that W^^ > everywhere (so long as the former is true, the latter can 

always be achieved by a suitable rescaling). One can therefore use the veto algorithm [T7] 
to correct down the emission rate. That is, whenever a trial has been selected according 
to eq. (E]), the ratio /W^^ in the chosen phase space point is the probability that this 
choice should be retained. If not, the evolution is continued downwards from the rejected 

scale {not from Q^ax)- This gives a change from eq. to 

^ZtZlQ^) = W^^iQ^) exp dQ- W^^iQ-)^ . (7) 

Note that, while all explicit dependence on W^^ is gone in eq. ([7]), an implicit dependence 
on the shower remains in two respects. Firstly, the Sudakov-factor modification of the 
basic ME shape reflects the order in which the shower algorithm sweeps over phase space, 
i.e. the shower definition. Secondly, if the as factors are omitted from the ly^^/iy^^ 
reweighting, the p5_"dependent expression used in the shower is retained, instead of the fixed 
value normal for ME's. 

Once the first emission has been considered, an uncorrected shower is allowed to continue 
downwards from the chosen scale. Thus, in this algorithm there is no fixed scale for the 
transition from ME to PS, but a smooth merging of one into the other. This cannot give 
rise to discontinuities in the behaviour at any phase space point, except of course at the 
soft/collinear shower cutoff. 

The formalism so far has only considered real emissions. For the e"'"e~ case, however, 
the cancellation of real (R) and virtual (V) divergences results in a particularly simple 
expression (when integrated over the possible orientations of the event) 

o-NLO = o-R + y do-R + (TV = (1 + ^ j cxb . (8) 
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It is therefore trivial to retain eq. ([7]) as an NLO prescription for the distribution of events, 
just by raising the cross section associated with each event from cjb to ctnlo- Identifying 
dQ"^ = dan/a-B, we arrive at the final equation for the differential cross section 

do- = o-NLO exp - / . (9) 

It is thus assumed that the 0{as) "new" part, o"nlo ~ <^b, of the total cross section 
should be associated with the same radiation function as ctb. Alternative choices on this 
count would only show up in 0{a'^), which is beyond the certified accuracy of the algorithm. 
Physicswise it is the minimal assumption, relative to having a different radiation function 
for the pure NLO part of the cross section. 

The merging formalism is easy to extend to the case of the emission of one extra gluon 
in any 1 — )■ 2 decay, and has also been worked out and implemented for all such cases 
within the Standard Model and the Minimal Supersymmetric Extension to the Standard 
Model [27j. It can also be apphed to ISR for 2 —t- 1 production of colour singlet particles, 
such as qq — )■ Z° [28]. Here, however, NLO corrections are not so trivial. Experimental 
practice has been to rescale the total cross section to the NLO answer, as in eq. thereby 
neglecting other potential kinematical differences between the LO and NLO answers. 

The merging approach has also been applied to HERWIG [29143T] . with two main dif- 
ferences. Firstly, since the HERWIG shower is ordered in angle rather than hardness, the 
ME/PS correction weight must be applied to any emission that is the hardest so far, rather 
than only to the first. Secondly, the HERWIG algorithm leaves holes in the phase space 
covered by the first shower emission, so it becomes necessary also to introduce a matching 
procedure, whereby such holes are filled directly by the ME rather than by the PS. 



2.2 The POWHEG strategy 

A general NLO cross section with hadronic incoming states will also contain remnant 
(counter) terms from the subtraction of initial-state coUinear singularities, which have 
already been incorporated into the PDF's. The complete differential cross section can 
therefore be written as the sum of contributions from leading order (Born), virtual, real 
and counter terms 

da = ddB + day + da^ — dac ■ (10) 

We note that the virtual term has the same n-body phase space as the Born term, while 
the real and counter terms have an [n + l)-body phase space. In the POWHEG method, 
the phase space kinematics are factorised in terms of Born (f ) and radiation (r) kinematic 
variables, such that the overall cross section may now be written as 

da = B{v) d*^ + V{y) d^^ + R{y, r) d$^ d$, - C{v, r) d$„ d$, . (11) 

Defining a function that integrates over the radiation variables 

B{y) = B{v) + V{y) + / d^r[R{y, r) - C{v, r)] , (12) 



once a Born event has been generated, distributed according to B{v) d$t,, the differential 
cross section for the hardest emission may be written as 



da = B{v)d^, 



B v) " \ Jv, B V 



(13) 
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The evolution variable of the POWHEG "shower" is taken as the kinematical p_i_ of the 
emission with respect to the parton that branches, identifying this highest p± branching as 
the hardest. For ISR this is the p± with respect to the beam axis. In this way, the hardest 
radiation is generated according to exact NLO matrix elements, but in a probabilistic, 
exclusive language, which can be directly interfaced to a suitable shower program. 

This expression shares many features with eq. The integral is from some upper 
scale associated with the Born event, with a lower cutoff to avoid soft/collinear regions. In 
both cases it is possible for an event to evolve down to this lower cutoff without radiating. 
Here, however, the constant NLO prefactor, ctnlo, is upgraded to B{v)d^y, a fully differ- 
ential quantity which will encapsulate all kinematical differences between the LO and NLO 
answers. It is thus in this term that all the sophistication and hard work of a full NLO 
calculation lies. As in the Pythia shower, the radiation is generated with a running as 
expression, evaluated at the p± scale of the emission. 

This formalism has been used for the hadronic production of vector bosons [32|l33] , heavy 
quark pairs [51], single tops [SS] and Higgs bosons via gluon/vector-boson fusion jSniEZ]- 
Codes for generating all these processes, except for vector boson pairs, is publicly available. 
The latest development is the POWHEG BOX; a general framework for implementing NLO 
calculations with the POWHEG method [38]. 

2.3 The Pythia transverse- momentum-ordered showers 

The Pythia 8.1 showers are ordered in transverse momentum [TB], both for ISR and for 
FSR. Also, multiparton interactions (MPI) are ordered in p± [39]. This allows a picture 
where MPI, ISR and FSR are interleaved in one common sequence of decreasing p± values 
[18] . This is most important for MPI and ISR, since they are in direct competition for 
momentum from the beams, while FSR (mainly) redistributes momenta between already 
kicked-out partons. 

The interleaving implies that there is one combined evolution equation 

dV _ f dVuPl ^ y^ dPlSR ^ y^ dpFSR A 

dp± \ dp± dp± dp± J 

that probabilistically determines what the next step will be. Here the ISR sum runs over 
all incoming partons, two per already produced MPI, the FSR sum runs over all outgoing 
partons, and p±max is the p_[_ of the previous step. Starting from a single hard interaction, 
eq. ( !T^ can be used repeatedly to construct a complete parton-level event of arbitrary 
complexity. Recently also rescattering has been included as a further (optional) component 
of the MPI framework [10] . 

The decreasing scale can be viewed as an evolution towards increasing resolution; 
given that the event has a particular structure when activity above some p± scale is resolved, 
how might that picture change when the resolution cutoff is reduced by some infinitesimal 
dpjj It does not have a simple interpretation in absolute time; all the MPI occur essentially 
simultaneously (in a simpleminded picture where the protons have been Lorentz contracted 
to pancakes), while ISR stretches backwards in time (and is handled by backwards evolution 
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and FSR forwards in time. The closest would be to view eq. ffn|) as an evolution 
towards increasing formation time. 

For the following, a relevant aspect is that the p± definition is not exactly the same for 
MPI, ISR and FSR in Pythia. For an MPI the p± is the expected one; the transverse 
momentum of the two scattered partons in a 2 —i- 2 process, defined in a Lorentz frame 
where the two incoming beams are back-to-back. To understand why problems arise for 
ISR or FSR, consider a q — > qg branching, where the p± of the emitted gluon is defined 
with respect to the direction of the initial quark. The p± as a function of the gluon emission 
angle 6 increases up till 90°, but then decreases again, — )■ for 6* — )■ 180°. Thus, an 
ordering in such a pj_ would classify a ~ 180° emission as collinear and occurring late in the 
evolution, although it would involve a more off-shell propagator than an emission at 90°. 
It could also erroneously associate a 1/pi divergence with the 6 — > 180° limit. Therefore it 
is natural to choose an evolution variable that does not turn over at 90°. 

To this end, consider a branching a ^ be (e.g. q — )■ qg), where z is defined as the 
lightcone (LC) momentum along the a axis that b takes. Then 

P±LC = ^(1 - ^Wa - (1 - ^Wb - zml , (15) 
and this equation can be used as inspiration to define evolution variables 

ISR : pi^^oi = (1 - z)Q'^ with ml = -Q^ and ml = ml = , (16) 
FSR : = z{l - z)Q^ with ml = and ml = ml = , (17) 

which are monotonous functions of the virtuality Q^. However, once a branching has been 
found and the kinematics is to be reconstructed, the interpretation is retained, but z is 
now replaced by a Lorentz invariant definition. For ISR this is chosen to be z = mf^/ml^, 
where r is the "recoiler", i.e. the incoming parton from the other side of the subcoUision. 
The actual p± of b and c then becomes 

Plb,c = (1 - ^)Q' pLvoI - , (18) 

where pl_evo\ max = (1 ~'2^)Qmax = i^^^Y^ar- Instead, for FSR, z = Eb/Ea in the rest frame 
of a and its colour-connected recoiler. The pl_^ ^ expression in this case becomes somewhat 
lengthier than eq. ( 1T8|) . but shares the same physical properties; at small values Pj_evoi 
pl_bc follow suit, and so correspond to identical l/pl_ singular behaviours, but the latter 
then turns around and vanishes when ]?5_evoi approaches the kinematical limit. 

2.4 Interfacing POWHEG-hvq to Pythia 

In this study, the interface from POWHEG-hvq to Pythia was performed using the Les 
Houches Event File (LHEF) [?T] format. Both programs come with a number of PDF sets 
available for selection and Pythia gives access to external PDF sets through the LHAPDF 
library [l2]. For consistency, the CTEQ6L [43j PDF set has been used throughout this 
study, as this set is available in all three programs used; no other differences in the PDF set 
itself have been taken into account. One caveat is the associated as running expression for 
a given PDF set. Pythia by default will use the first order expression, and allow as(^l) 
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to be set independently for spacelike and timelike showers. Here, as(M|) was set to match 
the one used in the PDF set, but with first order running. 

Historically the LHA conventions [44j, used by LHEF, only encompass one event "scale" , 
specified to be the factorisation one. The and Oem values at the renormalisation scale 
are also supplied, but not that scale itself, and there is no provision for shower-matching 
scales. A proposal for an extended standard that would store more scales has recently been 
presented |15]. When reading an LHEF into a shower program, the question then is: at 
what scale should the subsequent shower evolution begin? 

In Pythia, such choices need to be made for ISR and FSR separately. The basic 
principle should be to avoid doublecounting to the largest extent possible. With LHEF 
input it would be natural to let the showers begin at the factorisation scale, the only known 
one, and then proceed downwards. This is an allowed choice both for ISR and FSR, but 
for ISR another possibility is default. Here, events are split into two kinds, based on the 
absence or presence of particles that the shower can produce, i.e. d, u, s, c, b, g and 7. If 
the LHEF final state contains no such particles then the shower can populate the full phase 
space without any risk of doublecounting, which should give more realistic event shapes. 
This would be the case for many W/Z, top, Higgs and New Physics processes. If, on the 
other hand, the final state does contain particles that could be produced in the shower, then 
doublecounting would be more likely than not, and the factorisation scale again becomes 
the only reasonable choice. 

This more flexible attitude works well if the ME program does not mix different topolo- 
gies, but breaks down if, say, both tt and ttg/ttq events are supplied, with the latter 
intended to correspond to the fraction of tt events with an extra emission above the fac- 
torisation scale. When such a mixing is present, the showers off tt should not be allowed 
to populate the whole phase space. Thus, we conclude that interfacing to an LHEF cannot 
be done completely automatically, but must be made with some knowledge of which rules 
were used to produce the LHEF. 

Therefore, to provide a consistent interface to POWHEG, we must both consider its 
generation strategy and the information it stores in the LHEF. In events with an emission 
above the lower cutoff scale p±min, POWHEG chooses the factorisation scale to be the p± 
of the emission, p_lpowheg = P± (defined with respect to the beam axis in the case of 
POWHEG-hvq). 

When the POWHEG "shower" reaches p±min without any emissions, the factorisation 
scale is instead set equal to this p±min, P±powheg = P±min- Since p_Lmin normally is rather 
small, of the order of 1 GeV, the fraction of no-emission events is also small. 

The LHEF choice of storing and Oem rather than the renormalisation scale makes 
sense for traditional ME calculations, where typically one fixed scale is used. The POWHEG 
case is somewhat special, as it can have different renormalisation scales for the Born level 
process and for the subsequent radiative emission. 

For the case where POWHEG has already generated an emission, the obvious choice 
is to begin the evolution at the factorisation scale Pxpowheg; such that the shower will 
not generate harder emissions. However, there are then two potential complications: first, 
the mismatch between p_lpowheg and the lightcone-inspired p_Levoi scales of Pythia, and 
second, specifically for FSR, that p± is defined with respect to the direction of the emitting 
parton rather than to the beam axis. For ISR, eq. (|T8l) shows that pj_ < pxevoi, such that 
starting the shower from p_Levoi = P±powheg will lead to a small area of phase space not 
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Figure 1: Ratio of the kinematic p± of the first shower emission to the POWHEG emission, 
where the shower emission is (a) ISR or (b) FSR. In both cases, the results are shown 
when starting the shower at the factorisation scale and when starting the shower at the 
kinematical limit and vetoing above the POWHEG scale 

being covered. An FSR emission, on the other hand, may have a small p_Levoi with respect 
to the emitting parton but still a p± > p±powheg with respect to the beam axis. 

A simple solution to both these problems is instead to begin the shower at the largest 
possible scale, and then veto any emissions with a kinematic p± > p±powheg- If we consider 
the first shower emission, the multiplicative nature of the no-emission probability ensures 
that the emission rates below p±powheg will be correct, i.e. unaffected by the vetoes above 
it. The picture is slightly less clear for subsequent emissions; having accepted one shower 
emission below p±powheg, it is still possible for a later emission to be above it, since the 
first emission may well have had p_Levoi > P±powheg- The probability of such an occurrence 
is small, and effects formally of NNLO character, unenhanced by any large logarithms. 
They mainly show up for \ow-p± first emissions, where their importance on the event as a 
whole is less, but still nonzero. Another NNLO issue is that recoil effects from one emission 
can shift the p± of the previous ones, along with the hard process itself, either to lower or 
higher values. 

The current POWHEG-hvq generator uses a second order running as expression, but 
with a A fixed at rij = 5. Although slightly inconsistent, this only leads to changes beneath 
the Bottom and Charm scales. The A value is taken from a selected PDF set and is modified 
as in |l6]. In the LHEF output file, all incoming and radiated partons are massless and the 
values of the couplings, and aem, are set to zero in all events. 

To quantify how well the proposed interfacing works, we begin with top pair production 
(mt = 171 GeV), where all results are generated at LHC energies (pp, a/s = 14TeV). In this 
case, the number of light flavours, which defines the content of the proton and the allowed 
radiation flavours, goes up to and includes the bottom quark {rii = 5). To study the effect 
of the different shower starting scales, we examine the ratio of the flrst p± in the shower 
to the pj_ of the POWHEG emission (where the shower p± value is taken directly after 
the emission). This is shown in Fig. [H split into contributions from (a) ISR and (b) FSR. 
For ISR, we note that the ratios do not become larger than unity, but that when starting 
the shower at the factorisation scale, there is a region close to p±shower/p±hard = 1 where 
the phase space is not completely flUed. This gap is flUed when starting the shower at the 
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Figure 2: Final p± of the top pair for three different cases (see text). In the upper plot, it 
is difficult to distinguish the different curves, but the results are clearer to see in the lower 
plot, which shows the ratio to the "Kinematical Limit + Veto" result 



kinematical limit and vetoing emissions above p_lpowheg = P±hard; that is, the shower is 
running at "full steam" when it reaches the emission p± threshold. That the curve happens 
to be so flat near the endpoint is a coincidence, related to a cancellation between the blowup 
of the naive emission rate for smaller p_i_, with a Sudakov damping in the same limit, helped 
along by many events having a small pxhard in the first place. It does not happen for FSR, 
where the large top mass reduces radiation to a lower level overall. For FSR, when starting 
at the factorisation scale, there is a tail beyond unity, as discussed previously, while this no 
longer happens with the veto scheme in place. Note that the FSR rate below p±powheg 
does come up, however, since emitted partons now are allowed at a larger separation from 
their mother parton so long as they are still at small p± with respect to the beam. As 
already noted, the corrections are of higher order, but their inclusion is worthwhile for 
overall consistency. 

We move on to study the effect of shower emissions beyond the first. As discussed 
previously, after a first allowed shower emission below P±powheg; a subsequent emission 
may be generated above p±powheg- To examine this effect, three different cases are con- 
sidered: (1) the showers are started at the factorisation scale, (2) the showers are started 
at the kinematical limit and only the first shower emission is vetoed, and (3) the showers 
are started at the kinematical limit and all emissions above p±powheg are vetoed. As a 
measure of this effect, one may use the final p± of the top pair, which is shown in Fig. [2] 
for the three different cases. At first glance, all three approaches appear to give similar 
results, but the ratio plot reveals that the difference between the factorisation scale and 
the veto schemes is around 10%, while there is little difference between the two different 
veto schemes. Instead, in Fig. [3l we study the smearing of the first shower emission p± due 
to subsequent emissions. In this case only ISR is generated, so there is no ambiguity in 
picking out the first shower emission in the final event, and again, we show the ratio of the 
first shower emission to the POWHEG one, p±shower/p±hard, (a) immediately after the first 
emission and (b) after the shower evolution has finished. Fig. [3^ shows the same features 
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Figure 3: Ratio of the kinematic p± of the first shower emission to the POWHEG emission 
for ISR only, (a) Shows the resuhs immediately after the first emission, while (b) shows 
the results after the full shower evolution. In both cases, the "Kinematical Limit + 1st 
Emission Veto" and "Kinematical Limit + Veto" curves lie on top of each other 
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Figure 4: Ratio of the kinematic p± of the first shower emission to the POWHEG emission 
for bottom pair production, where the shower emission is (a) ISR or (b) FSR 

as Fig. [H although with different normalisation due to the lack of FSR. In Fig. [3]d, the first 
shower emission p± is smeared by small amounts due to subsequent emissions, but again 
the difference between vetoing just the first emission and all emissions is negligible. 

For bottom production (mt = 4.8 GeV), the number of light flavours is now set to four 
{ill = 4), meaning that there will be no incoming or radiated b quarks. This is different 
from the default Pythia settings, where b quarks are both taken from the beam and 
allowed to be created in radiative emissions. The results are shown in Fig. IH again split 
into contributions from (a) ISR and (b) FSR. Although the overall pattern of radiation is 
different, due to the smaller bottom mass, the details with respect to the interface with 
Pythia show the same features as Fig. [TJ 

2.5 Poor Man's POWHEG 

The POWHEG approach is a very powerful tool for NLO studies. Unfortunately, 
POWHEG-based implementations are not available for all relevant processes. Many older 
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NLO calculations are instead available in terms of phase space slicing results (see e.g. |17]). 
That is, the (n + l)-body cross section is based on the pure real-emission cross section da-R, 
without any Sudakov corrections. Unresolved emissions, virtual corrections and countert- 
erms are lumped together with the Born cross section to provide the n-body cross section. 
The absence of a Sudakov factor means that soft /collinear divergences in do"R are not damp- 
ened, thus requiring a reasonably large cutoff scale to stay away from the region where cr^ 
would turn negative. On the other hand, duR should be used where it can. Therefore an 
and (Tn+i typically are chosen to be of the same order, while cr„ <C cr„+i in POWHEG, 
where a much smaller cutoff can be used without any risk of inconsistencies. 

In order to make these older NLO codes also useful to the experimental community they 
need to be interfaced to event generators, for further showering, MPI and hadronisation. It 
is then convenient to add Sudakovs to dan, to bring it closer to the POWHEG approach. 
This would also smoothen out the transition between (n + 1)- and n-body phase space, 
in the spirit of modern-day leading-order matching procedures, already mentioned in the 
introduction. 

As a practical example we mention the Wgamma_nlo package for W7 production to 
NLO ^S], with up to one additional quark or gluon in the final state, and allowing for 
anomalous WW7 couplings. Here an implementation of the principles to be presented is 
already available j49] . 

Assume a lowest-order process that does not contain any coloured particles in the final 
state. The NLO processes thus include one quark or gluon extra, with divergences when 
its p_L — > 0. A scale pxmin is used to separate the (n -|- 1)- and ra-body phase space. An 
{n + l)-body state is characterised by the transverse momentum of the quark or gluon, p^^, 
with > p_Lmin- We now want to include the Sudakov to express that is the hardest 
jet of the event, i.e. that there are no jets at a larger scale. 

In the CKKW-L [2] approach a fictitious trial shower is used to provide this. The better 
this shower attaches to the correct ME behaviour, the more accurate the Sudakov will be 
(thus, an advantage to having a more accurate default behaviour of the shower, as we strive 
for in this article). The point should not be overstressed, however; at large p± the Sudakov 
suppression is negligible anyway, and at small p± the universal behaviour should dominate. 

The {n + 1) topology has to be projected down to the core n process, without the 
extra parton, to provide the starting point for the trial shower. This requires a choice 
of which side the emission occurs on. Sometimes flavours allow only one possibility; for 
ug — )• W"'"7d the hard subprocess must be ud — )■ W+7 since the gluon cannot couple to 
the W'^7 state. For ud — )■ W"'"7g the emission could be on either side, and the relative 
probability for a shower emission is related to the respective 1/pievoi- Since z = m^/m^j^^ 
is independent of emission side, it follows that Pj_evoi oc oc 1 =F cos^, where 9 is the angle 
of the gluon in the rest frame of the subcoUision. The constancy of z also implies that 
the splitting kernel values are the same on both sides, and so give no net contribution. In 
another case, where the core process would be different depending on the side of emission, 
the weight of the respective splitting kernel would have to be taken into account, but still 
at the same z. In summary, the relative probability for a gluon emission from side 1 is 
Qi'^/iQi'^ + Q2^) = Q\I{Q\ + Qi) = (l + cos^)/2. The choice of side determines a new 
x' = xz for it, while x on the other side is unchanged. A combined rotation and boost can 
bring the n particles to this new frame. 

From there, the fictitious shower is allowed to generate the first/hardest emission at a 
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scale . The probability for pj^ > is precisely the desired Sudakov, in the CKKW-L 
spirit. That is 

(i) if > the original n + 1-body topology is retained; 

(ii) else the projected n-body topology is selected. 

In case (i), a normal shower can then be started up from p^^ and downwards; that is, 
further jets may be generated above the p±min scale, to give (n + 2)-body topologies etc. 
The events in case (ii) are to be lumped together with the ones already originally classified 
as n-body, and allowed to shower from the p_Lmin scale downwards. 

In this approach, there is not guaranteed to be a smooth matching at p^min, at least 
for the p±_ scale of the hardest emission. The hope is that this step will be smeared out by 
subsequent showers and hadronization. Also note that, unlike POWHEG defined by eq. f|T3|) . 
but in line with MC@NLO, the high-p^ tail is defined by the LO (n + l) expression, without 
any B{v)/B{v) "K factor". 



3 POWHEG without NLO 

3.1 Correction to the Pythia shower 

As will be shown in Section 1X51 for top pair production, an evolution with Q^^x — ^ (power) 
overestimates the high-pj_ tail while Q^ax = "^It (wimpy) underestimates it, cf. |2T]. This 
is not so surprising, as follows. Let us recall that the matrix element for a QCD process 
such as gg — )■ ggg scattering roughly behaves like [501 l5T] 



da ~ ^(2) (p^, + + p^3) , (19) 

pii p±2 pis 

where the p±i are the transverse momenta of the three outgoing gluons. In the limit 
P±3 P±iy where |p±i| ~ |p±2|5 this reduces to l/pliPia, i-e. the hard interaction behaves 
like dp5_i/pli and the subsequent shower emission of an additional gluon like dp^^/pi^. 
Put another way, for a fixed p±i, you may distinguish a low-p_L3 region with a fall-off like 
dpis/pis and a high-p_|_3 region where the fall-off instead is like dpi^/p^^, with a smooth 
transition when p±3 ^ p±i. In practice you would not want to simulate the process like this, 
of course, but reserve the hardest propagator to be described by matrix elements, which 
also would avoid doublecounting problems. 

Obviously the picture is not equivalent for the gg — )■ ttg, but let us apply a similar 
reasoning to the extra gluon in this case. For small p±g the picture of an ISR branching 
g — )■ gg followed by a hard process gg — )■ tt ought to be a vahd approximation, and so we 
expect a dp^g/p5_g falloff. At large p±g, on the other hand, it would make more sense to 
think in terms of an ISR branching g — )■ tt followed by a hard process gt — )■ gt, and thus 
a shape more like dp5_g/plg- Now, we don't simulate the latter kind of hard processes, e.g. 
because the top is so heavy that a top PDF is not a fruitful approximation, and so we 
would like to obtain these high-p^g configurations in the context of the simulation of the 
gg — )• tt process. This leads to an ansatz of the form 

dPisR ±_ k'M' 

dpi pi km^ + pI ' ^ ' 
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where is a reasonable scale to associate with the hard 2—7-2 process (discussed further 
below) and k"^ is a fudge factor of order unity, parameterising at what scale the transition 
from a to a l/p± behaviour occurs. 

How generic would such an ansatz be? First of all, for QCD processes involving light 
quarks and gluons, the showers should be cut off at (around) the scale of the hard process, 
or else one would doublecount, since in this case all possible 2 — )• 2 hard subprocesses one 
could construct out of a 2 — )■ 3 process are already being simulated. Secondly, for 2 — )■ 1 
production of a colour singlet particle, the dp\/p']_ ansatz works well up to the kinematical 
limit |28]. A reasonable assumption is that this generalises to processes where two or more 
colour singlet particles are produced in the core process, while a shape similar to eq. (120]) 
should occur in processes that involve one or several coloured particles in the final state. 

The argument for such a difference is one of colour coherence; with colour charge in 
both the initial and the final state one expects a destructive interference between ISR and 
FSR emissions that limits the radiation [52] , while no such interference occurs with colours 
only in the initial state. One can then argue exactly what the scale appearing in 
eq. ( 12U]) should be; for pair production of coloured particles, the factorisation or renormali- 
sation scale should be a reasonable choice, but when considering the production of a mixed 
coloured/non-coloured final state, by the coherence argument, it is primarily the coloured 
particles that should play a role in this scale. We note that the default choice for internal 
Pythia 2 — 2 processes is to set the renormalisation scale equal to the geometric mean 
of the squared transverse masses of the two outgoing particles, while the factorisation scale 
is set to the smaller of the squared transverse masses of the two outgoing particles. In the 
next sections we will check how well the ansatz of eq. ( 120|) fares for a number of different 
processes. 



3.2.1 Hard emission probability 

In order to study processes for which an implementation of POWHEG is not yet available, 
we turn to the MadEvent matrix element generator. In eq. f[T^ we see that the POWHEG 
Sudakov, used to generate hard emissions, contains only Born and real terms; all NLO 
corrections are contained in a separate prefactor. 

We therefore use MadEvent to generate events with an extra jet in the final state in order 
to extract the cross section for real jet emission, dcr/j. An approximate "POWHEG" style 
probability for emission is then formed by normalising to the overall lowest order (Born) 
cross section (generated by simulating the corresponding 2 — )■ 2 process in MadEvent) and 
adding a Sudakov. This leads to a probability distribution described by 



where the p± integral begins at the kinematical limit. As noted previously, the NLO prefac- 
tor can lead to kinematical differences in the resulting distribution. Under the assumption 
that such kinematical differences are small (to be addressed further in the case of top pair 
production), a qualitative comparison to the Pythia shower can be made. The effect of 
the Sudakov will be most visible in the low-p^ region, where the ME da^ calculation is 



3.2 MadEvent 




(21) 
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divergent. Here the integrand will blow up and the Sudakov will make the distribution 
turn over, such that it will have a unit integral (up to cutoff effects). Given that the low- 
p± region is where the parton shower should be most accurate, we expect the turnover 
in the Sudakov-modified MadEvent distribution to roughly correspond with the Pythia 
distribution. 

All Standard Model masses in MadEvent are set equal to the default Pythia values. 
Variations in particle widths between MadEvent and Pythia are not expected to play a 
large role and are therefore neglected. The definition of both the proton and jets includes 
bottom quarks, as would be included in the Pythia default settings. In what follows, we 
will only consider heavy final states (top mass and higher), such that we can neglect FSR 
and only consider ISR. In this case, there are only soft and coUinear singularities from the 
branching of incoming partons. A \ow-p_i_ cutoff of 2GeV is introduced on real-jet emission, 
such that these divergences are avoided. No further cuts are applied. 



3.2.2 Scale corrections 

In generating events with MadEvent, there are again some different choices available relating 
to running ctg expressions and PDF's. The choice of PDF determines the running of ctg and 
the value as{Mz) and as before, the CTEQ6L PDF set was chosen. The ME calculations 
were generated with fixed renormalisation and factorisation scales, taken to be the geometric 
mean of the masses of the two heavy final-state particles. The results of these choices remain 
evident in the ME-derived jet distribution. In eqs. ([2]) and (|3]) it was shown that the Pythia 
shower algorithm instead picks scales related to the the evolution variable when generating 
emissions. These differences can be non-negligible, especially in the low-px region, where ag 
will become large and we expect a similar turnover to the Pythia distribution. To account 
for these effects, additional weights are applied to MadEvent distributions after generation, 
with the new scale taken to be the kinematic p± of the jet. 

A correction for is simple to achieve, weighting events by the ratio of the new and 
old tts values. With M defined as the fixed renormalisation/factorisation scale used in the 
ME calculations, we have a weight 

^ML . (22) 

With ISR, correcting for differences in the factorisation scale is more difficult. The distribu- 
tion for real emission from the ME calculation contains two PDF factors from the incoming 
partons, unlike ISR generated from the shower algorithm. This difference makes it unclear 
what correction should be applied, especially as in our ME calculation, we do not know 
which incoming parton has branched. Taking an exclusive point of view, there should have 
been no emissions between and p\ on either side of the event, and so both should be 
reweighted, overall giving a correction factor 

Xlfl{Xi,pl) X2f2{x2,pl) . . 

Xi/i(Xi,M2) X2/2(X2,M2) ' ^ ^ 

This may be a slight over-correction, however, and one would expect the results with just 
the tts factor and those with this additional PDF factor to bracket the "correct" distribution. 
We examine this further in Sec. 13.31 for the case of top pair production. 



15 




(a) (b) 

Figure 5: Two possible graphs for squark pair production with an additional jet which 
do not directly correspond to a 2 —t- 2 process and a Pythia shower splitting 

3.2.3 Graph topologies 

When moving to a final state containing an extra jet, MadEvent will correctly generate all 
possible topologies, but when comparing against Pythia, some of these graphs may not 
correspond to a shower history. Two examples are given in Fig. |5] for the case of up squark 
pair production in the MSSM. 

In (a), the two t-channel squark propagators cannot directly be reproduced by a 2 — )■ 2 
process and a shower emission. In a traditional merging approach, these topologies can 
always be associated to some 2 — )■ 2 hard process with the radiation associated to an 
incoming leg. As radiation off a massive leg is suppressed, topologies of this kind are not 
expected to dominate the overall behaviour of the real-emission matrix elements. 

Instead, in (b), the squark and jet are created through the resonant decay of a gluino. 
This resonant mechanism leads to events with a different topology than we wish to con- 
sider; as Pythia works in the narrow width approximation, where the gluino production 
and decay are described as part of the separate 2 — )• 2 process ug — )■ lig, it would be 
doublecounting to include the same graph as a 2 — )■ 3 process. 

To remove double counting of this type, we use the method outlined in [33]; the idea 
here is to assign propagators to events statistically, based on the relative size of the relevant 
squared amplitudes. Such events then explicitly contain this propagator information in the 
LHEF output, meaning they can be easily vetoed when subsequently being processed. As 
a simple check that the results are reasonable, many processes were also generated with 
these resonant graphs manually removed, e.g. as in [21]. Although this procedure is not 
fully gauge invariant, as the widths of the resonances are small, we expect the results to be 
in line with the MadEvent veto method. 

3.3 Standard Model processes 

We begin by comparing the jet emission probability in top pair production of POWHEG, 
MadEvent and the first Pythia shower emission (ISR only). In all studies that follow, 
events are generated at LHC energies (pp, y/s = 14 TeV) and, additionally, when the 
damping ansatz is used, the scale of eq. fl20|) is set to be the factorisation scale unless 
stated otherwise. For top production, we restrict both Pythia and MadEvent to QCD 
production only, as this is the dominant contribution to the cross section, and also allows a 
consistent comparison to POWHEG. As noted previously, the default Pythia ISR shower 
will begin evolution at the kinematical limit, given that a top cannot be produced in the 
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Figure 6: First emission probability as a function of p± for top pair production in POWHEG 
and MadEvent. Ratio plots are normalised to the POWHEG result 



shower. The POWHEG to MadEvent comparison allows us to check the validity of the 
MadEvent approximation method outlined in Sec. 13.2.11 The comparison to Pythia then 
lets us examine the effectiveness of the ansatz described in Sec. 13.11 

The probability distributions for jet emission in top pair production are expected to be 
roughly in agreement over the entire p± range for POWHEG and MadEvent, as long as the 
kinematical differences that come from the full NLO prefactor of eq. f|T3|) are small. This 
is shown in Fig. El separately for low and high p±_ regions, for three different sets of the 
MadEvent results. All three have the Sudakov correction applied, but one set additionally 
has the as scale correction ("os corrected"), and the other both the as and factorisation 
scale corrections ("as & /i corrected"). These will all meet when the p± scale matches that 
of the fixed factorisation/renormalisation scale used in the MadEvent generation, but begin 
to diverge as the p± rises and falls away from this value. In the low-p^ region, the Sudakovs 
make the distributions turn over, while in the high p±_ tail, there is good agreement between 
the distributions of POWHEG and the approximation from MadEvent; in the ratio plot 
the POWHEG curve sits between the two different corrected MadEvent curves, but overall 
closer to the as & /i-corrected one. 

We now move on to compare POWHEG against the Pythia shower and the ansatz 
of Sec. 13.11 The expectation is that the default Pythia power shower will not fall off 
quickly enough and therefore overestimate the emission probability in the high p± tail. In 
Fig. [7] we show, again separately for low and high p± regions, the POWHEG results against 
four variations of the Pythia shower (wimpy, power, damped with k = 1 and damped 
with k = 2). Again, there is good agreement in the low-p^ region. In the high-p^ tail, 
the wimpy shower sits clearly below the POWHEG result, while the Pythia power shower 
does indeed overestimate the jet emission probability. The damping procedure, particularly 
with k = 2, then brings the Pythia distribution into closer agreement with POWHEG. 

It is interesting to study the kinematical differences that come from the full NLO pref- 
actor in POWHEG by considering the top pair rapidity both with and without such a 
prefactor. When full NLO corrections are included, we not only expect an overall change 
in cross section, but also that there may be kinematic differences in the distributions. In 
Fig. [8^, the ratio of two POWHEG pair rapidity distributions is shown. The first is from 
the default POWHEG-hvq generator, with a full NLO prefactor, while the second is gen- 
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Figure 7: First emission probability as a function of p± for top pair production in POWHEG 
and Pythia, split into (a) the low-p^ region where the Sudakov makes the distributions 
turn over and (b) the high-p_L tail 




Figure 8: (a) Ratio of top pair rapidity between POWHEG with the full NLO prefactor, B, 
of eq. ([13]) and B only. Note the suppressed zero on the y axis, (b) Jet rapidity distributions 
of Pythia and MadEvent for p±j^^ > lOOGeV, normahsed to the POWHEG result 



erated by replacing the B term of eq. f[T^ with just the Born contribution B. The figure 
shows both features; an overall shift in cross section (a fc-factor of around 1.5) and a small 
shift in pair rapidity to more central regions. 

We finally compare the jet rapidity distributions from POWHEG, MadEvent and 
Pythia. In Pythia, the rapidity of the extra jet is taken immediately after the first 
shower emission (if present). For MadEvent, however, we can not apply the Sudakov cor- 
rection in this case, and the direct output of the generator is taken. We restrict ourselves 
to p_Ljet > lOOGeV, such that the effect of the Sudakov and other scale corrections will be 
small. Here it is harder to know what to expect. For the damped Pythia shower, it is 
clear that by only reducing the high-p_|_ tail relative to the power shower, events will be 
taken out of the central region. These distributions are shown in Fig. [Sb, normalised to 
the POWHEG result. MadEvent and POWHEG are in good agreement, although with an 
overall shift in cross section. The default Pythia shower is strongly peaked in the central 
region, but comes into better agreement when damped. 
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Figure 9: First emission probability as a function of p± for Z pair production. Results show 
Pythia compared to the approximate MadEvent {as & /i corrected) prescription 

(a) (b) 
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Figure 10: W pair production, (a) Shows the invariant mass of the jet and W boson when 
the jet is a bottom quark for three different MadEvent runs (see text). The large resonant 
peak in the "MadEvent" sample has been cut off above 0.15. (b) Shows the first emission 
probability as a function of p± 

We now study processes where there is not (yet) a "correct" POWHEG result to compare 
with, but instead rely on the MadEvent approximation. If the ansatz of Sec. 13. H is correct, 
the damping of the Pythia shower is not expected to help improve the p_i_ distribution in Z 
and W pair production in a significant way. Here, there are no coloured final-state particles 
to use as a guide in setting the scale of eq. (|2U]) . so we continue to use the factorisation 
scale when generating the damped results. Again, in this case, the Pythia default is to 
start ISR shower at the kinematical limit. Fig. [9] shows the results for Z pair production 
for the default, damped [k = 2) and wimpy shower against MadEvent {as Sz /x corrected). 
The results show that the default shower, although giving a slightly too hard p± tail, still 
does a reasonable job of reproducing the MadEvent curve without additional damping. 

The case of W pair production is useful as a check on the MadEvent resonance veto 
method (see Sec. I3.2.3p . When bottom quarks are allowed in both the incoming beams 
and in the definition of the extra jet, there is now a resonant contribution coming from 
t — >■ bW. In Fig. [TOk . the invariant mass of the jet and the matching W are shown when 
the jet is a bottom quark. A clear peak is visible when those events with a resonance are 
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Particle 


Mass (GeV) 




561.1 


g 


607.7 


X? 

xf 


96.7 
181.7 



Table 1: Subset of MSSM SPSla masses used in process generation 

not vetoed ("MadEvent"), while this peak is almost completely removed when these events 
are vetoed ("MadEvent + Veto"); there is only a small remaining trace of it from e.g. 
interference terms. As an additional check, also plotted is the result when resonant graphs 
are explicitly removed from the MadEvent generation ("MadEvent + Graph Removal"). 
Fig. [TOb . shows the results for the first emission probability in the p± tail for the MadEvent 
vetoed sample. As a final check, the radiation pattern was compared to a sample not 
including bottom quarks, and found to contain no large differences. 

3.4 MSSM studies 

To further study the production of heavy final states, MSSM processes were chosen split into 
three groups: coloured final states, non-coloured final states, and coloured/non-coloured fi- 
nal states. All events were generated with the SPSla [51] set of parameters (the relevant 
final-state masses used are shown in Tab. [1]). Here, also the lowest order 2 — )■ 2 processes 
were generated using MadEvent, and fed into Pythia for showering, again with the fac- 
torisation and renormalisation scales fixed at the geometric mean of the masses of the two 
heavy final-state particles. Note that this is slightly different from the default Pythia 
internal 2 — 2 process scale choices and that this choice enters into the scale selection 
of the damping ansatz, eq. [201 

3.4.1 Coloured final states 

For the fully coloured final states, all events were generated with QCD only; as with top 
pair production, this is where the dominant contribution to the cross section lies. As before, 
the expectation is that the coloured final states will benefit from a damping of the high-p^ 
tail. Figure fTTl shows the tail of the p± distributions, all compared to the MadEvent as & 
fi corrected data, for (a) ulUl, (b) Uig and (c) gg. Similarly to top pair production, the 
power shower overestimates the high-p_L tail, while the damping ansatz brings the curves 
into better agreement. For ulUl production, k = 2 leads to the best agreement, while for 
ul g and g g production, the MadEvent curve lies between the damped k = 1 and k = 2 
curves. 

3.4.2 Non-coloured final states 

We move on to study the non-coloured final states, Xi Xi aiid Xi Xiy where we expect the 
damping ansatz not to improve the pj_ tail of the parton shower. Here, for both processes, at 
lowest order, there are large resonant H'^/A° contributions. These are relatively long-lived 
intermediate states which we expect to follow the rules for 2 — )■ 1 singlet production, where 
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Figure 11: First emission probability as a function of pj_ for (a) ulUl, (b) Uig and (c) gg 
production. In (a) the k = 2 and MadEvent curves lie on top of each other 




Figure 12: First emission probability as a function of p± for (a) Xi Xi and (b) Xi Xi 
production 



the shower already does a good job in covering the entire phase space. For both processes 
then, the MadEvent veto scheme is used to remove events of this type. 

The results are shown in Fig. [12] for (a) Xi Xi (b) Xi Xi production. The results 
follow the pattern for Z/W^ pair production; in both cases, the default power shower does 
a reasonable job in the high-p^ tail. 
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Figure 13: Pythia first emission probability for tlie damped sliower {k = 2) in UlXi 
production, wliere is set eitlier to the factorisation scale or to the mass of the heavy 
outgoing squark. The default power shower curve is shown for comparison 

3.4.3 Coloured/non-coloured final states 

Finally, we study mixed coloured/non-coloured final states. Here, we would expect some 
measure of damping to improve the shower description due to the presence of colour in 
the final state. As discussed in Sec. 13. one issue here is the meaning of the scale of 
eq. (120|) . With the mixed final state, the difference between the factorisation scale and e.g. 
just the mass of the final-state coloured object is large enough to give noticeable differences 
in the damped shower tail. As an example. Fig. [T3] shows the difference in the damped 
shower [k = 2) in ulXi production, when is set to the factorisation scale and when it 
is instead set to the mass of the outgoing squark. In the remaining results of this section, 
we retain the choice of setting to the mass of the outgoing coloured state. Also note 
that all of the processes studied in this section have strong resonant contributions, where 
a squark decays into a neutralino/chargino and a jet, making it difficult to generate large 
statistics when using the MadEvent veto scheme. 

Results are shown in Fig. [H] for (a) ulXi (b) ulXi- While the curve for ulXi 
does show the expected behaviour, with the damping improving the shower description, for 
UlXi "^^ obtain the expected results. Instead, Fig. \T5k shows the results for vllXi 

when resonant graphs are manually set to zero in the MadEvent generation code. Here, 
beyond p_L ~ 1 TeV, the tail does appear to fall off with a damped behaviour. Comparing 
the real-emission Feynman diagrams for niXi against those of ulXi ^ a key difference is the 
appearance of right-handed intermediate squarks. In Fig. [T5b . we show the results when 
the right-handed squark masses {ur, d^, and cr) are set high. In this case, we do recover 
the behaviour of the ulXi result. 

Finally, in Fig. [16], we study the processes (a) gx? and (b) gXi- Again the expected 
behaviour is not apparent; both sit closest to the default power shower curve. 

In summary, for the mixed processes studied in this section, our ansatz is not an obvious 
improvement relative to the power shower. There is a nontrivial dependence of the emission 
pattern on the SUSY parameter choices that we do not understand, and do not go on to 
study further at this time. 
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Figure 15: First emission probability as a function of for ulX?, where in (a) resonant 
graphs have been removed and (b) the right handed squark masses have been set high 




Figure 16: First emission probabihty as a function of p± for (a) gXi and (b) gXi 
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4 Summary and Outlook 



One direct outcome of this study is an improved matching between the existing POWHEG 
program family and the Pythia p^-ordered showers. The main point is that the Pythia 
p± definition does not quite agree with the normal kinematical one, although the two do 
agree closely over most of the phase space populated by showers. In hght of this, the 
numerical effect of these improvements are fairly modest, but, in view of the important role 
that we foresee for POWHEG-based studies at the LHC, there is every reason to have the 
interface to Pythia well understood. 

The POWHEG approach, in complete or in simplified form, also allows us to test the de- 
fault behaviour of the p^-ordered showers. In particular we address the issue of "power" vs. 
"wimpy" showers, i.e. what starting scale to use for the downwards evolution, specifically 
for pair production of heavy particles. Here we show that many processes obey an inter- 
mediate behaviour, where the characteristic dp5_/pi shower fall-off is replaced by a steeper 
dpx/pi fall-off for large p± values. That is, emissions are allowed up to the kinematical 
limit, but at a dampened rate. The damping can be approximated by a factor 

^'dampen - ^2^2 + p2 ' (24) 

for the emission of a parton of transverse momentum p±, where is a characteristic scale 
of the hard process and is a fudge factor. 

More specifically, we have seen that for coloured pair production, with set equal to 
the factorisation scale (noting the slightly different scale choice used in the Pythia internal 
processes and the external MSSM ones) and = 2, we get a reasonable shower behaviour 
for the different final states studied here. Instead, when the particles in the final state of 
the hard process are colour singlets only, we have seen that there is no need to impose 
a damping, consistent with previous results for single W/Z production. The differences 
between these two cases can be understood based on the destructive interference between 
ISR and FSR that is expected for coloured particles but not for the uncoloured ones. In no 
case is the wimpy shower a good choice. 

For the mixed coloured/non-coloured final states of Sec. 13.4.31 we have argued that 
should instead be related to the coloured partons only, due to this coherence argument. 
Unfortunately, for these processes, the results are not as expected. We do not currently un- 
derstand the structure of the real-emission matrix elements that gives rise to the behaviour 
seen here and, specifically, whether it is an artifact of the complicated SUSY structure or 
an inherent property of the radiation pattern. 

Not studied here is the case of the production of light coloured particles, i.e. of normal 
QCD jets. There, the hard process and the showers produce the same kind of partons, and 
doublecounting becomes a main concern. We intend to return to this class of event. 

With the advances in computational tools, allowing automatised (Born-level) higher- 
order calculations, one may question the need for more accurate showers. We believe there 
are two main points in favour of improving showers, especially when this can be done with 
only a modest effort. One is that new physics scenarios are continuously being proposed, 
where higher-order calculations may be overkill for first studies, but nevertheless a realistic 
population over all possible event topologies is useful, even if off by a factor of two or so 
in some tails. The other is that trial showers as a means of obtaining Sudakov factors 



24 



is a crucial ingredient of CKKW-L matching schemes, so that the quahty of matching is 
improved if the quahty of the shower is also improved. 
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